Dynamic generation of spin orbit coupling 
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Spin-orbit coupling plays an important role in determining the properties of solids, and is crucial 
for spintronics device applications. Conventional spin-orbit coupling arises microscopically from 
relativistic effects described by the Dirac equation, and is described as a single particle band effect. 
In this work, we propose a new mechanism in which spin-orbit coupling can be generated dynamically 
in strongly correlated, non-relativistic systems as the result of fermi surface instabilities in higher 
angular momentum channels. Various spin-orbit couplings can emerge in these new phases, and 
their magnitudes can be continuously tuned by temperature or other quantum parameters. 

PACS numbers: 71.10. Ay, 71.10. Ca, 71.10. Hf 
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Most microscopic interactions in condensed matter 
physics can be accurately described by non-relativistic 
physics. However, spin-orbit (SP) coupling is a notable 
exception, which arises from the relativistic Dirac equa- 
tion of the electrons pj. The emerging science of spin- 
tronics makes crucial use of the S-P coupling to manip- 
ulate electron spins by purely electric means. The pro- 
posed Datta-Das device Q modulates the current flow 
through the spin procession caused by the SP coupling. 
More recently, Murakami, Nagaosa and Zhangjj, 3 pro- 
posed a method of generating the dissipationless spin cur- 
rent by applying an electric field in the p-doped semicon- 
ductors. This effect and the similar proposal for the n- 
doped semiconductors [5| both make crucial use of the SP 
coupling. In contrast to the generation of the spin current 
by coupling to the ferromagnetic moment, purely elec- 
tric manipulation has an intrinsic advantage. However, 
unlike the ferromagnetic moment, which can be sponta- 
neously generated through the strong correlation of spins, 
the conventional wisdom states that the SP coupling is a 
non-interacting one-body effect, whose microscopic mag- 
nitude is fixed by the underlying relativistic physics. 

On the other hand, recent interests are revived on the 
Landau-Pomeranchuk (L-P) [(| fermi surface instabili- 
ties, largely in connection with electronic liquid crystal 
states with spontaneously broken rotational symmetry 
0- 0- S lit! ilTL IT^ I , and in connections with hidden or- 
ders in heavy fermion systems [l3L IT3. ITH l . Varma's recent 
work showed that the L-P instability could lead to the 
opening of an anisotropic gap at the fermi surface |l3|. 
In this paper, we show that the SP coupling can be gen- 
erated dynamically in a non-relativistic system through 
strong correlation effects as the L-P instability in the 
spin channel with higher orbital angular momentum. It 
emerges collectively after a phase transition, which is con- 
tinuously tunable either by temperature or by a quantum 
parameter at zero temperature. Unlike the ferromagnet, 
our ordered phase keeps time reversal symmetry. Also in 
contrast to the L-P instabilities considered by the major- 
ity of previous theories, most translationally invariant liq- 
uid phases in our model do not break rotational symme- 
try, and some of them preserve time reversal and parity 
symmetries as well. Most correlated phases in condensed 



matter physics are characterized by their broken sym- 
metries |16). Solids break translational symmetry, liq- 
uid crystals break rotational symmetry, superfluids and 
superconductors break gauge symmetry and ferromag- 
nets break time reversal symmetry and rotational sym- 
metry. As far as we are aware, the new phase reported in 
this work is the only one besides the fermi liquid which 
does not break any of the above symmetries. It is distin- 
guished from the fermi liquid by only breaking the "rel- 
ative spin orbit symmetry" , a concept first introduced in 
the context of the 3 He liquid |17| . 

We first discuss the dynamic generation of SP coupling 
from the L-P instability within the Landau-Fermi liquid 
theory triggered by the negative Landau parameter f 
and then present its exact definition. This instability 
lies in particle-hole channel with total spin one and rela- 
tive orbital angular momentum one. Operators in matrix 
forms are defined as Q^ a {r) — ^Pai^a^^—iW^tp^r), 
where Greek indices denote the direction in the spin 
space, Latin indices denote the direction in the orbital 
space, and the operation of V a on the plane wave is de- 
fined as V<V^ = (V7|V|)e* F = fcV'K Q» a (r) is 
essentially the spin-current operator up to a constant fac- 
tor. We use a Hamiltonian similar to that of Ref. but 
in the channel: 

H = Jd 3 r Vi(r)(e(V)- M )^(r) + h^Q^lf) 

+ \ f d 3 rd 3 ? f?(r-r)Q^(r)Q» a (r<), (I) 

where \i is the chemical potential and the small h^ a is 
dubbed as the "spin-orbit field" , which plays a role sim- 
ilar to the external magnetic field. For later convenience 
, we keep both the linear and the cubic terms in the ex- 
pansion of the single particle dispersion relation around 
the fermi wavevector kf, e(k) = VfAk[l + b(Ak/kf) 2 ], 
with Afc = k — kf. We assume that the Fourier com- 
ponents of fi(r) take the form /f (q) = J dr e lqr /f (r) = 
fi /(l + K\f{ L \q 2 ) and define the dimensionless Landau pa- 
rameter F® — Nff®, where Nf is the density of states at 
fermi energy. The symmetry of the Hamiltonian Q is a 
direct product SO(3)l <8> 50(3)5 m the orbital and spin 
channels. 
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We define the spin-orbit susceptibility as Xna,vb = 
{Qiia)/h v b in the limit — » 0, which is diagonal, i.e. 

= X$iJ.a,5 v b, in the normal fermi liquid phase. The 
fermi liquid correction to \ is given by 

m* 1 

Xfl — Xo— * . [,„„ ■ (2) 

TO 1 + -T-f /O 

with to* /to the ratio between the effective and bare 
masses. The spin-orbital susceptibility is enhanced for 
Fi < and divergent as the critical point F" — —3 is 
reached. 

In the mean-field (MF) analysis, the p-h channel triplet 
order parameter is defined as n M ° (f) = — J df /" (r — 
r')(Q fla (r')) ) and the external spin-orbit field h^ a is set to 
zero. Using the uniform ansatz ra M<l (r) = n Ma , Eq. ^is de- 
coupled into H M f = J (Pr^{r){e(V) - ^V^-iV") - 
At) V> (f ) + Vn^ a n» a / ( 2 1 /f | ) , with V t he space volume . The 
self-consistent equation for the order parameters reads 

/d 3 k 
" ^KW)), (3) 

which is valid when the interaction range tq = y «|/f | 
is much larger than the distance between particles 1/kf, 
i.e., the dimensionless parameter rj — kfy/K\ff\ 1. 

The phase structures can be determined from the 
Ginzburg-Landau (G-L) free energy, which is similar to 
the trip let pairing order parameter in the 3 He system 
[T7L Hq . Under the independent SO(3) rotations in the 
orbital and spin spaces Rl and Rs, n^ a transforms as 
7i^ a — > RL,^n y ' b R^\ a . Furthermore, is even under 
the time-reversal but odd under the parity transforma- 
tion. With these symmetry requirements, the G-L free 
energy can be constructed up to the quartic order as 

F(n) = A tr[n T n] + B 1 {tv[n T n]) 2 + B 2 ti[(n T n) 2 ]. (4) 

Compared with the complex order parameter in the su- 
pcrfluid 3 He case, the reality of the n pa restricts the free 
energy to contain only two quartic terms. Explicitly, 
tr[n T n] = n^nP* = X 2 + \\ + A|, and tr[(n T n) 2 ] = 
n^ a n^ h n va n vh = \\ + A| + A|, where A 2 2)3 are eigen- 
values of n T n. For B 2 < or B 2 > 0, Eq. 0] favors 
the the structures of (A?, A 2 ., A 2 ) to be proportional to 
(1,0,0) or (1,1,1), respectively. We name them as a 
or /3-phases whose general order parameter matrix struc- 
tures are given by 

n na _ f n d^e a a-phase, for B 2 < , g s 
\ n F>^ a /3-phase, for B 2 > 0, 

where d and e are two arbitrary unit vectors in the spin 
and orbital space respectively, D^ a is any SO (3) rotation 
matrix, and n is a real number. In other words, the cor- 
relation functions of operators Q^ a acquire a long range 
part in the ordered states 

<q-(^v)> - >< { X: (6) 




a-phase (3-phase 

FIG. 1: The L-P instability in the F£ channel, with dashed 
hnes marking the fermi surface before symmetry breaking. 
The fermi surface distortion is anisotropic in the a-phase, 
while it is isotropic in the /3-phase with dynamic generation 
of spin-orbit coupling. 

as \r — r*| — » oo. This correlation function gives the 
rigorous definition for the new phases, independent of 
the approximate Fermi liquid theory used here. 

The a-phase is a straightforward generalization of the 
nematic fermi liquid $j to the triplet channel as shown 
in Fig. ^ where the spin and orbital degrees of free- 
dom remain decoupled, and the rotational symmetry 
is broken. Taking a special case n^ a = n8^ z 5 a z, the 
dispersion relations for spin up and down branches are 
£, A {k)i, 2 = e(k) — fi =F ncos6> respectively, where 9 is the 
angle between k and z-axis. The fermi surfaces for the 
two spin components are distorted in an opposite way as 
Ak fh2 (6)/k f = ±xcos8 (1 - bx 2 cos 2 9) - 1/3 x 2 , with 
the dimensionless parameter x — n/(vfkf). The chem- 
ical potential /i is shifted to ensure the particle number 
conservation as 5/j,/(vfkf) — —x 2 /3. The remaining sym- 
metry is 50(2)^ <g> SO(2)s with the Goldstone manifold 
5£ (g) Sg. Two Goldstone modes are the oscillations of 
the distorted fermi surfaces, and the other two are the 
oscillations of the spin directions. 

In the /3-phase, the rotational symmetry is preserved 
with the dynamic generation of spin-orbit coupling as 
shown in Fig. ^ F° r example, with the ansatz n Ma = 
n8fj, a , the MF Hamiltonian is reduced to 

H MF =^ t (fc)(e(fc)-//-n5 ! -fe)V>(fc). ( 7 ) 

k 

The single particle states can be classified according 
to the eigenvalues ±1 of the hclicity operator a ■ k, 
with dispersion relations £ B (fc)i,2 = — /x ± n. 

The fermi surface distortions of two helicity bands are 
Afe/i^/fc/ = ±x(l — bx 2 ) — x 2 and the chemical poten- 
tial shift Sn/(vfkf) = —x 2 . Similar to the superfluid 
3 He-B phase, the /3-phase is essentially isotropic. The 
orbital L and the spin 5 angular momenta are no longer 
separately conserved, but the total angular momentum 
J = L + 5 = remains conserved. The Goldstone mani- 
fold is [SO(3) L ®SO(3)s]/SO(3) L+s = 50(3) with three 
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branches of Goldstone modes. For the general case of 
n fj,a = nDua j it is equivalent to a redefinition of spin op- 
erators as S'^ = S v D va 8 a ^, thus fermi surface distortions 

remain isotropic and J' = L + 5" is conserved. A similar 
symmetry breaking pattern also appears in the quantum 
chromodynamics where the two-flavor chiral symmetry 
SU(2) L x SU(2) R is broken into the diagonal SU(2) L+R 
[lfij . In that case, both SU(2)l,r are internal symme- 
tries, and thus there is no flavor-orbit coupling. 

The coefficients of the G-L free energy Eq. 0] can be 
microscopically derived from the MF theory as 
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where b describes the cubic part of the dispersion e(k), as 
explained earlier. With the definition of 8 = 1/|-Ff| — 1/3, 
the L-P instability takes place at 8 < i.e., Ff < —3. 
For b < 1/3, i.e. B2 < 0, the a-phase appears with 
1*12 _ 



For b > 1/3, i.e. B 2 > 0, the /3-phase 

— The largest fermi 



l"i _ 2(B 1 +B 2 ) 

is stabilized at with Inl 2 



2(3Si+S 2 ) 

surface distortion Akf tTnax /kf in the a-phase is larger 
than the uniform one Akf jkf in the /3-phase, thus a large 
positive b is helpful to the /3-phase. However, we empha- 
size that this is only one of the options to change the sign 
offla. 

To apply the a and (3 phases in the lattice system, we 
only need replace the 50(3)l symmetry with the lattice 
point group. For example, for the simple cubic lattice, 
we define Q^ a = i{c> (x)^ 'c{x + e a ) — h.c.} with e a the 
base vector in the a-direction. The unbroken symmetry 
is Ol ® SO(3)s where Ol is the orbit lattice octahedral 
group. The mean field Hamiltonian for the /3 phase reads 



Hmf = ^2 ^!( £fc ~ V~ na » sin k a8^a)^k, 



(9) 



with lattice momentum k restricted in the first Brillouin 
zone. The helicity structure for each k is aligned along 
the direction of (sin k x , sin k y , sin k z ) , which breaks the 
symmetry down to the combined octahedral rotation in 
the orbit and spin space Ol+s- As a real space analogy, 
the hexagonal non-colinear anti-ferromagnet YMnC>3 [2£j, 
EH ] has the spin order pattern inside the unit cell which is 
also invariant under the combined spin-orbit point group 
rotations. The difference is that the spin order in the j3 
phase lies in the momentum space and no spatial spin 
order exists each lattice site. The lattice a-phase was 
originally studied under the name of "spin-split state' by 
Hirsch 23 to explain the phase transition at TV = 311-ftT 
in the Chromium system. 

By reducing the space dimension to 2, the mean 
field Hamiltonian for the /3-phase reduce to the famil- 
iar Rashba[2^ and Dressselhaus Hamiltonians [24j in the 
2D semiconductor heterostructures. The order parame- 
ter rifxa is a 3 x 2 matrix with fi = x,y, z and a — x, y. 



Its third row of fj, = z can be transformed to zero by 
performing suitable SO (3) rotation on the index /x, thus 
we take n Ma as a 2 x 2 matrix. The G-L free energy is 
also the same as in Eq. QJ but with the new coefficients 



2 k |/f| 2 ) 



B\ = 



N, 



32 u jk j 



B 2 = 



bN f 
Svjkj 



,(10) 



and the L-P instability occurs at F£ < —2. The a and /3- 
phase structures are similar as before in Eq. [S] However, 
there are two options in the /3-phase with n^ a — nD^a, 
where D^a is a 0(2) matrix. If det£> = 1, then J z — 
L z + S z is conserved. With the MF ansatz n 
we arrive at the Rashba-like Hamiltonian 



fi.a 



ne 



H R = d 



>t| e (v)-ne^ a ^HV°)}V- (11) 



If det D = — 1 , J z is not conserved while the en- 
ergy spectrum and fermi surface distortions are still the 
same as the case of det D — 1 . With the MF ansatz 
n^a = n diag{l,— 1}, we arrive at the 2D Dresselhaus- 
likc Hamiltonian as 



d 2 r 



V> f {e(V) - n[a x (-iV x ) - a y (-iV v )] }^>. (12) 



If we generalize the mechanism of dynamical genera- 
tion of the spin-orbit coupling to the spin 3/2 fcrmionic 
system, an interesting phase can be obtained which pre- 
serves all familiar symmetries including the parity sym- 
metry, breaking only the relative spin orbit symmetry. It 
has been recently shown that any generic model of spin 
3/2 with local interactions has an exact SO(5) symme- 
try in the spin space [25|. The four spin components form 
the spinor representation of the 5*0(5) group. Using the 
Dirac T matrix defined there, the spin 3/2 Landau in- 
teraction functions are classified into the SO(5)'s scalar, 
vector and tensor channels p5| : 

f af 3,ys(p,p) = fipj) + f v (pJ)(T a /2) aP (r a /2) 75 

+ /*(p, P ')(r Q V2)a /3 (r a V2) 75 . (13) 

We further pick out its L = 2 part of the orbital an- 
gular momentum in the spin 2 vector channel denoted 
as the F% channel. We define operators Q fia (r) — 
i>l(r)Kpd a {V)Mr) (1 < M.a < 5), where d Q (V) = 
[VSV x V y , -VSV X V Z , V3V y V z , -±(3V 2 -V 2 ), ^(V 2 - 
V 2 ,)]. The model Hamiltonian is constructed as follows: 

H= f d 3 r ii(f)(e(V)-n)^p a (f) 

+ If dVdV tt{?-r*)Q» a {r)Q» a (r>), (14) 

with the symmetry of SO(3) L <E> SO(5) s . The order 
parameter is defined as before n A ' a (r) = — J df^f^r — 
r'){Q^ a {r')) and the L-P instability occurs when F% = 
iV//|(g = 0)<-5. 
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The ordered phases after the L-P instability can also 
be classified into two categories as before: the a-phase 
with anisotropic fermi surface distortions and /3-phase 
with spin-orbit coupling. The detail phase structures 
are much more complicated here. For example, the 
a-phase has two non-equivalent configurations because 
the L = 2 channel fermi surface distortions can be ei- 
ther uniaxial or biaxial. A comprehensive classification 
of all the possible phases is quite involved and is de- 
ferred to a future work. We focus here on the /3-phase 
with the order parameter structure n^ 0, — fi8^ a . In 
this case, the MF Hamiltonian is reduced into Hmf — 

J d 3 r^ (r){e(V) - \i- nS fia T^d a ('V)\ip(r). From the re- 
lation between the matrices and the quadratic form 
of spin 3/2 matrices 5 p|, it can be easily recognized the 
Luttinger-like Hamiltonian |26| 

H L = J cPripi (r){e(V) -fi- fi(-iV • S) 2 }tp(r), (15) 

which is the standard model for the hole-doped III- 
V semiconductors. The original symmetry SO(3)l <8> 
SO(5)s in Eq. HH is broken into SO(3) L +s with ten 
branches of Goldstone modes. Different from the spin 
1/2 case, the order parameter n^ a here is even both un- 
der the time-reversal and parity transformations. 

It is remarkable that most standard spin-orbit sys- 
tems, including the familiar Rashba, Dresselhaus, and 
Luttinger-like Hamiltonians can all be generated dynam- 
ically by our mechanism, which can compliment the con- 
ventional mechanism based on relativistic single parti- 
cle physics. If a system has a small spin-orbit coupling 
to begin with, it acts like a spin-orbit field h va which 
rounds off the phase transition discussed here. One way 
to experimentally detect the dynamic generation of the 



spin-orbit coupling is through the anomalous Hall effect 
(AHE), which relies on the spin-orbit coupling. There- 
fore, detecting the AHE signal turning on at a phase 
transition directly demonstrates that the phase transition 
breaks the relative spin-orbit symmetry. We propose to 
systematically search for these new phases in 3 He, ultra- 
cold atomic systems, semiconductors, heavy fermion ma- 
terials and ruthenates, both in experiments and in mi- 
croscopic numerical simulations. The new /? phases can 
be experimentally detected through the standard predic- 
tions based on the spin-orbit coupling, but would have 
the remarkable and distinctive feature that these effects 
turn on and off near a phase transition. We conjecture 
that besides the familiar superfluid A and B-phases, 3 He 
may contain the new phases proposed in this work. The 
Landau parameter F£ in 3 He was determined to be neg- 
ative from various experiments [27l l28l I29I l3p| such as 
the normal-state spin diffusion constant, spin- wave spec- 
trum, and the temperature dependence of the specific 
heat etc. It varies from around —0.5 to —1.2 with in- 
creasing pressures to the melting point, reasonably close 
to the instability point F" = —3. Even though we pre- 
sented mean field descriptions of the new phases with 
dynamically generated spin orbit coupling, the existence 
of these phase can obviously be studied by exact micro- 
scopic calculations of the correlation function Eq. © for 
realistic models. 
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